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Abstract 

Exact solution of Schrodinger equation for the pseudoharmonic potential is ob- 
tained for an arbitrary angular momentum. The energy eigenvalues and correspond- 
ing eigenfunctions are calculated by Nikiforov-Uvarov method. Wavefunctions are 
expressed in terms of Jacobi polynomials. The energy eigenvalues are calculated nu- 
merically for some values of I and n with n < 5 for some diatomic molecules. 
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1 Introduction 



The study of anharmonic oscillators in three dimensions has much attention in chemical 
physics phenomena [1-45]. They have wide applications in molecular physics [1]. The solu- 
tions of Schrodinger equation for any £-state for such potentials have also a special interest. 
The Morse potential is commonly used for an anharmonic oscillator. However, its wave- 
function does not vanish at the origin, but the Mie type and the pseudoharmonic potentials 
do. The Mie type potential has the general features of the true interaction energy [1], in- 
teratomic and inter- molecular and dynamical properties of solids [2] . The pseudoharmonic 
potential may be used for the energy spectrum of linear and non-linear systems [3]. The 
Mie type and pseudoharmonic are two exactly solvable potentials other than the Coulombic 
and harmonic oscillator. This potential may be considered as an intermediate potential be- 
tween the harmonic oscillator potential and anharmonic potentials, as the Morse oscillator 
(MO) potential, which is a more realistic potential, in good agreement with the experimen- 
tal spectroscopical data [47]. A comparative analysis of potentials HO-3D (3-dimensional 
harmonic oscillator potential) and PHO is performed in [2]. 

To solve the Schrodinger equation for the pseudoharmonic potential, the Nikiforov- 
Uvarov method [46] which is introduced for the solution of the hypergeometric type second 
order differential equations appeared in the time-independent problems. The method is 
based on the determination of the solution in terms of special orthogonal functions for any 
general second-order differential equations [47]. 

By taking an appropriate s = s(r) coordinate transformation, the Schrodinger equation 
in spherical coordinates is reduced to a generalized equation of hypergeometric type as 

where, a(s) and a(s) are polynomials, at most second-degree, and f(s) is in a first-degree 
polynomial. 

The contents of this work is as follows. In section 2, solution of Schrodinger equation 
for any £— state is introduced. In section 3, solutions for Mie-type and pseudoharmonic 
potentials are given. As an example, some numerical values of the energy levels of N 2 , CO, 
NO and NH molecules are computed. Section 4 is devoted for the calculations. 



2 Method and Solutions 

The pseudoharmonic potential has the form [3] 

where Vo is the dissociation energy and tq is the equilibrium intermolecular separation. The 
radial part of the Schrodinger equation for a diatomic molecule potential is 

/ h 2 1 d ( 2 d \ £(£ + l)h 2 \ „ „ , x 
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where p is the reduced mass of the diatomic molecules, n denotes the radial quantum num- 
ber (n and t are named as the vibration- rotation quantum numbers in molecular chemistry), 
r is the internuclear separation. Substituting the explicit form of V(r), we get 



#JW " + ^ + % (** - * f" - ^ - C ^Pf\ B-M - 0. (4) 



o?r 2 r dr % \ \r r ) 2pr 

By defining a coordinate transformation 

r 2 = s, (5) 
and the following dimensionless parameters 



we obtain 



2 Vo P 

e = (E ne + 2V ) ^ (7) 



Thus, comparing the Eqs. (1) and (8), we obtain the following polynomials 
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f(s) — —, a(s) — s, and a(s) = —a 2 s 2 + es — (3. (10) 

To find a particular solution of Eq. (9), we apply separation variables as R n i{s) = <f>(s)y(s), 
Schrodinger equation, Eq. (1) is reduced to an equation of hypergeometric type 

a(s)y" + T(s)y' + Xy = 0, (11) 

where <f>(s) is a solution of (f>'(s) ^ <f)(s) = ir(s)/a(s) and y(s) is a hypergeometric type 
function whose polynomial solutions are given by the well-known Rodrigues relation 



y(s) = y n (s) = [<j n {s) P {s)\ . (12) 

Here, a(s) is a polynomial of degree at most one, p(s) is a weight function and B n is a 
normalization constant. It satisfy the condition 

{ap)' = rp. (13) 
The function n and parameter A are defined as 
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and 



(14) 



\ = k + n'. (15) 

Here 7r(s) is a polynomial with the parameter s. Determination of k is the essential point 
in the calculation of n(s). fc is determined by setting the square root as the square of a 
polynomial. Thus a new eigenvalue equation for the Schrodinger equation becomes 

X = X n = -nr'- n{n ~ 1) a"; n = 0,l,2..., (16) 

where 

r(s) = f(s) + 2tt(s), (17) 

and it will have a negative derivative. Therefore the polynomial of ir(s) is found in four 
positive values 



= -±±(as-y/p + ±), k_ = e- 1) (18) 



where 7 = 2a^Jp+ J. By following Eq. (17), we get 



, -J ± (as + , /or £;+ = e + 7 

r( S ) = ^ + 2 (19) 
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"I ± (« s - >JP + 5) , for k_=e- 7. 



To have a negative derivative of r(s) and a physical eigenfunction, we take the following 
substitution: for fc_ = £ — 7, we have 



t(s) = -7" \ as -\ P + T?\i (2°) 



4 \ V 16 

and 



r(s) = l-2(as-J/3+4l, (21) 
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Thus, using Eq.(16) we obtain 



A = A n = 2na, (22) 
and using Eqs.(22) and (15) we obtain the following energy eigenvalues, 
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= I 2n • 1 • 2, >' • — | n. 



(23) 



Substituting Eq.(23) into Eq. (7), the energy eigenvalue becomes 



h 2V 

E n i — — 2Vo H \ 

r V /i 



(2n + l) + 2, 



+ 16 



(24) 



In order to calculate the wave function, the weight function p can be obtained by using 
Eq.(13) as 



p(s) = s 2 ^+¥e e -2- (25) 
and using Eq.(12) which is Rodrigues relation, we can get the function y, 

y n ( s ) = s i+2 V / ^¥ e ~ 2as 
Finally the wave function is obtained as 



(26) 



(27) 



3 Conclusions 

We have studied analytical solution of Schrodinger equation for a Pseudoharmonic potential. 
Energy eigenvalues and the corresponding wave functions are calculated for a diatomic 
system with any angular momentum L . The NU method is used in the computations. 
Numerical values of energy for A^2, CO, NO and NH molecules are calculated for different 
values of principal and angular quantum numbers n and i. Results are listed in Table I. 
The potential parameters are obtained from Ref. [10]. 
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Table 1: Energy eigenvalues (in eV) of pseudoharmonic potential for N2, CO, NO and CH 
diatomic molecules with different values of n and £. Necessary parameters are given in [10]. 



n 


t 


iV 2 


CO 


NO 


CH 








0.10915590 


0.10193061 


0.08248827 


0.16863440 


1 





0.32734304 


0.30567217 


0.24735916 


0.50500718 




1 


0.32784167 


0.30615078 


0.24778171 


0.50859034 


2 





0.54553018 


0.50941373 


0.41223005 


0.84137996 




1 


0.54602881 


0.50989234 


0.41265260 


0.84496312 




2 


0.54702603 


0.51084953 


0.41349768 


0.85212458 


4 





0.98190446 


0.91689685 


0.74197183 


1.51412550 




1 


0.98240309 


0.91737546 


0.74239438 


1.51770870 




2 


0.98340031 


0.91833265 


0.74323946 


1.52487010 




3 


0.98489606 


0.91976835 


0.74450700 


1.53560020 




4 


0.98689026 


0.92168247 


0.74619689 


1.54988430 


5 





1.20009160 


1.12063840 


0.90684272 


1.85049830 




1 


1.20059020 


1.12111700 


0.90726527 


1.85408150 




2 


1.20158750 


1.12207420 


0.90811035 


1.86124290 




3 


1.20308320 


1.12350990 


0.90937789 


1.87197290 




4 


1.20507740 


1.12542400 


0.91106778 


1.88625710 




5 


1.20756990 


1.12781650 


0.91317990 


1.90407610 
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